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Abstract. Far-infrared maps obtained with ISOPHOT have been searched for 
point-like sources. The majority of the 55 sources is believed to be extragalactic 
and in most cases no previously known sources can be associated with them. Based 
on the far-infrared spectral energy distributions it is likely that dust-enshrouded, 
distant galaxies form a significant fraction of the sources. 

We estimate the number densities of extragalactic sources at 90u.m, 150u.m 
and 180um wavelengths and at flux density levels down to ~100mJy. The counts 
are compared with models of galaxy evolution. The counts exceed the predictions 
of current models, even those with strong evolution, and no-evolution models are 
rejected at a high confidence level. 

Comparison with recent results from COBE mission indicates that at 90u.m 
the detected sources correspond to <;20% of the extragalactic background light. At 
longer wavelengths the corresponding fraction is 10%. 



1 Introduction 

In the far-infrared the cosmic infrared background (CIRB) consists of radi- 
ation emitted by galaxies, intergalactic gas and dust, photon-photon inter- 
actions (7-ray vs. CMB) and, possibly, by decaying relic particles. A large 
fraction of the energy released in the universe since the recombination epoch 
is contained in the CIRB. An important aspect is the balance between the 
UV-optical and the infrared backgrounds. With recent observations at in- 
frared and sub-mm wavelengths it has become obvious that star formation 
efficiencies in the early universe derived from optical and UV observations 
are underestimated (e.g. Madau et al. jl7); Steidel et al. |2^| ). 

The analysis of the data from the COBE DIRBE (Hauser et al. §) and 
FIRAS (Fixsen et al. indicated a FIR CIRB flux at a surprisingly high 
level of ~1 MJy sr _1 between 100 and 240 urn. Similar results had been 
obtained already by Puget et al. |pj and Schlegel et al. |^2). Because of 
the great importance of the FIR CIRB for cosmology these results require 
confirmation by independent measurements. 

The final goal of the ISOPHOT CIRB project is the determination of the 
FIR CIRB flux level. First steps are the measurement of the CIRB fluctua- 
tions and the detection of the bright end of FIR point source population. The 
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Table 1. The studied fields. The columns are: name of the field, coordinates of 
the centre of each field, area of the map, name of the ISOPHOT filter used in the 
observations, number of raster positions observed, step between adjacent raster po- 
sitions in the staring mode mapping and integration time. The distance of adjacent 
scans was in all cases identical to the raster step used along the scan line 



Field 


Map Centre 
RA(2000.0) DEC(2000.0) 


A 

Area 
(sq.degr.J 


Filter 


Rasters 


Step 
(arcsec) 


tint 

(s) 


VCN 


15 15 21.7 


+56 28 58 


0.030 


C_90 


10x4 


90 


46 










C.135 


10x4 


90 


46 










C.180 


10x4 


90 


46 


VCS 


15 15 53.1 


+56 19 30 


0.023 


C.90 


21x2 


90 


46 










C.135 


21x2 


90 


46 










C.180 


21x2 


90 


46 


NGPN 


13 43 53.0 


+40 11 35 


0.27 


C.90 


32x4 


180 


23 










C.135 


32x4 


180 


27 










C.180 


32x4 


180 


27 




13 42 32.0 


+40 29 06 


0.53 


C.180 


15x15 


180 


32 


NGPS 


13 49 43.7 


+39 07 30 


0.27 


C.90 


32x4 


180 


23 










C.135 


32x4 


180 


27 










C.180 


32x4 


180 


27 


EBL22 


02 26 34.5 


-25 53 43 


0.19 


C.90 


32x3 


180 


23 










C.135 


32x3 


180 


27 










C.180 


32x3 


180 


27 


EBL26 


01 18 14.5 


01 56 40 


0.27 


C.90 


32x4 


180 


23 










C.135 


32x4 


180 


23 










C.180 


32x4 


180 


23 



ISOPHOT CIRB project has potential advantages over the DIRBE analysis: 
(1) with the much smaller f.o.v. ISOPHOT is capable of looking at the darkest 
spots between the cirrus clouds; (2) in spite of the smaller f.o.v. ISOPHOT's 
sensitivity surpasses that of DIRBE in the important FIR window at 120 - 
200 |xm; (3) with the good spatial and spectral sampling the galactic cirrus 
can be separated. 

We have mapped four low-cirrus regions at high galactic latitude at the 
wavelengths of 90, 150, and 180 urn. Here we report on the point sources 
(galaxies) found in the FIR maps. The FIR source counts are important for 
the study of the star formation history of the universe and for the testing of 
models of galaxy evolution. 

The point source extraction is based on the fitting of the detector footprint 
to spatial data. The method is different from those used in most previous 
studies (e.g. Kawara et al. @; Puget et al. pofl ) where the source detection 
has been based on the analysis of the detector signal as a function of time. Our 
analysis is therefore independent of and complementary to previous results. 
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2 Observations 

The observations were performed with the ISOPHOT (Lemke at el. Jl3[) 
aboard ISO (Kessler et al. |flj). The maps were made in the PHT22 staring 
raster map mode (see Table]]]). The area covered is ~1.5 square degrees. The 
fields have low surface brightness and in some cases there is some redundancy 
i.e. the observed pixel rasters partly overlap each other. 

The data were processed with PIA (PHT Interactive Analysis) versions 
7.1 and 7.2. The flux density calibration was made using the FCS (Fine 
Calibration Source) before and after each map. The accuracyof the absolute 
calibration is expected to be better than 30% (Klaas et al. fll2|j). 

Data reduction from the ERD (Edited Raw Data) to SCP (Signal per 
Chopper Plateau) was performed using the pairwise method (Stickel p5|). 
Instead of making linear fits to the ramps consisting of the detector read- 
outs one examines the distribution of the differences between consecutive 
read-outs. The mode of the distribution is estimated with myriad technique 
(Kalluri & Arce |^]) and is used as the final signal for each sky position. The 
pairwise method is robust against glitches and the analysis was based on the 
data reduced with this method. 

3 Detection procedure 

The source detection was performed in two steps using data processed to 
the AAP (Astrophysical Applications Data) level with PIA and the pairwise 
method. The data consists of surface brightness values with error estimates. 
Flat fielding was performed with custom routines. 

Each surface brightness value was compared with the mean of the region 
within a radius of ^three times the size of the detector pixel. Values more 
than 0.7a above the local background were considered as potential point 
sources. A model consisting of a point source and a constant background was 
fitted into each region surrounding a candidate position. Footprint matrices 
were used to calculate the contribution of the point source to the observed 
surface brightness values. The free parameters of the fit were the source flux 
density, the two coordinates of the source position, and the background sur- 
face brightness. The formal errors are used to calculate the probability, P, 
that the detection is not caused by background noise. 

The completeness of the source detection and the number of false de- 
tections were studied with simulations. The results were used to adjust the 
probability level that was used for discarding uncertain detections. The ratio 
p between the source flux density and the background rms noise, o^g, was 
used as a criterion to discard uncertain sources. Since p is not directly re- 
lated to the probability obtained from the footprint fit it can be used as an 
additional safeguard against false detections. 
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The actual source list consists of sources detected at two or three wave- 
lengths with a spatial distance between detections of less than 80". The ex- 
pected number of co-incidental associations is no more than ~10% and does 
not significantly affect the source counts. 

4 Source counts 

The surface density of sources is estimated by dividing the number of detected 
sources at a given flux density level with the corresponding 'effective' map 
area. 

As the first approximation the effective area corresponding to a given 
flux density level was taken to be the sum of those maps where sources with 
equal or lower flux densities were detected. The cumulative source densities 
obtained at 90uxn, 150u.m and 180u.ni are shown as histograms in Fig. |[ 
Two sets of sources were used. The first set consists of all detections (dotted 
line) while in the second set there are only sources detected at more than one 
wavelength (solid line). The results are similar for 150u.m and 180u.m and 
significant differences are seen only at 90 urn. This is expected, since sources 
seen at f50u.ni are likely to be seen also at 180um (and vice versa) while 
more of the 90um sources remain unconfirmed at the longer wavelengths. 

A third estimate for the cumulative source densities (dashed lines in 
Fig. ||) was obtained by selecting sources based on the ratio p (see Sect. |[). 
For each flux density level sources with p > po were selected and no confir- 
mation at other wavelengths was required. The effective area was obtained 
by integrating the total area where the rms noise was below l//?o times the 
source flux density. The values of po were selected based on simulations. At 
the bright end the results agree with earlier histograms since no sources are 
rejected and the corresponding areas converge towards the total area mapped. 

At very low flux density levels the small number of sources leads to large 
uncertainties. The values obtained below f50mJy for 150u.m and f80u.ni are 
probably only indicative. 

5 Discussion 

5.1 Cirrus confusion 

We have checked the probability that some of the sources detected are small 
scale cirrus structures (cirrus knots). Results of Herbstmeier et al. show 
that at the scale of the C100 beam size, d ~ 45", the expected cirrus fluc- 
tuation amplitude is below 10 Jysr" 1 / 2 for all our 90u.m maps. According to 
Gautier et al. [Q this corresponds to a flux density of 4 mJy which is clearly 
below the flux densities of the faintest 90 urn detections. Therefore, cirrus 
is not likely to be a significant contaminant in the source counts although, 
because of the non-gaussian nature of the cirrus fluctuations (Gautier et al. 
[Q), our source list may still contain a few cirrus knots. 
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5.2 Comparison with galaxy spectra 



In Fig. [l] we compare the average of the source spectra with the spectra 
of the galaxies Arpl93 and NGC4418. In the sample of luminous infrared 
galaxies presented by Lisenfeld, Isaak & Hills |u| Arp 193 has the lowest and 
NGC4418 the highest estimated dust temperature. 

In the rest frame the SEDs of luminous infrared galaxies reach maxima 
between 60/im and 100/xm (Silva et al. [p3l ; Devriendt, Guiderdoni & Sadat 
[0; Lisenfeld, Isaak & Hills (IB)). The spectra of our FIR sources are flat 
in the observed wavelength range and the emission peak is typically above 
90^im. This is consistent with most sources being at redshifts 0.5^2^1. 

The FIR emission maximum of normal spiral galaxies is also located close 
to 100/im (e.g. Silva et al. 1 23 ) . However, in the case of a spiral galaxy a 
FIR detection at the level of 0.1 Jy would correspond to a visual magnitude 
brighter than 16 and the optical counterpart should be readily visible. The 
lack of clear visual counterparts indicates that most of our sources are likely 
to be more distant and luminous infrared galaxies. 
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Fig. 1. The average of the source spectra and the model spectra of Lisenfeld et 
al. ([^|) for the luminous infrared galaxies Arp 193 and NGC4418. The dashed 
lines show the spectra of the two galaxies shifted to z=l.Q. The flux density scale 
is arbitrary 



5.3 Comparison with previous source counts 

From Fig. || we obtain lower limits of 9xl0 4 sr _1 , 2.0 x 10 5 sr _1 and 2.2 x 10 5 sr _1 
for the cumulative source counts >100mJy at 90pm, 150 pm and 180 pm, 
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respectively. Similarly, above 200 mJy the source densities are 1.3xl0 4 sr 1 , 
6xl0 4 sr" 1 and 7xl0 4 sr~ 1 . The dotted histo gram, which lies mostly between 
the two other curves, gives at 100 mJy values 1.4x10 sr , 2.5xl0 5 sr~ 1 and 
S^XlC^sr -1 at 90 urn, 150 urn and 180 nm. 

Kawara et al. 10 have performed similar counts using ISOPHOT observa- 
tions of the Lockman Hole. The cumulative source counts were l.lxl0 5 sr _1 
and 1.3x10 s sr -1 at 95 (i.m and 175 \im, respectively, for sources 5,y>150mJy. 
At 90 urn our results are close to their values. Our lowest estimates that were 
based on detections at two or three wavelengths are, however, lower by more 
than a factor of two. This could indicate that the requirement of having de- 
tections at more than one wavelength does indeed underestimate the counts. 

Latest results from the ELAIS survey (Oliver et al. fl8| ), covering an 
area of 11.6 square degrees at 90|^m, are similar to the preliminary results 
reported by Rowan-Robinson et al. f2l| . At ^150mJy level the source density 
is 3.3x 10 4 sr -1 i.e. a factor of three lower than Kawara et al. |ici[| . The result 
is close to our lowest estimate as shown in Fig. 0. 

Linden- V0rnle et al. 14 have performed source counts in a 0.4 square 
degree field at 60|J.m and 90|xm. At 90p.m the results are at 100 mJy level 
similar to the ELAIS counts but are lower at other flux density levels. Also 
in their re-analysis of the Lockman Hole field they found significantly lower 
source densities than reported by Kawara et al. |l0| . 

Puget et al |Q have published source counts based on 175|J.m ISOPHOT 
observations of the Marano field covering ~0.25 square degrees. The source 
densities (e.g. ~4-10 5 sr _1 at 120 mJy) are slightly higher than our counts. 

The source counts obtained in these different ISOPHOT studies are mainly 
within a factor of ~2 from each other but direct comparison is made difficult 
by the differences in the source flux calibration adopted by the different au- 
thors. We have found a ~30% difference between the DIRBE surface bright- 
ness values and our data calibrated with PIA version 7.3 (Juvela et al. B). 
This uncertainty is indicated in Fig. 3 by the horizontal lines that are used 
to represent our source counts. At 90|J.m the DIRBE calibration results in 
lower flux densities and at 150|J.m and 180|J.m in higher flux densities. These 
values are better for the comparison with e.g. the results of Oliver et al. jl8| 
at 90|J.m since they used DIRBE as the basis of their calibration. 



5.4 Comparison with galaxy models 

The source counts at 150um and 180|xm are much higher than predicted by 
no-evolution models (e.g. Franceschini et al. At 180 |J.m the difference is 
a factor of five and these models can be safely rejected. Our results indicate 
that the luminosity or the number of galaxies must evolve strongly with z. 

Guiderdoni et al. || have presented semi-analytic models for the galaxy 
evolution. In their model E (which predicts the highest source counts) both 
the 'burst' mode of star formation rate and the relative number of ultralumi- 
nous IR galaxies (ULIRGs) increase with z; at z = 5 half of all galaxies are 
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ULIRGs. The model is in good agreement with extragalactic background light 
measurements in both optical and infrared. The model predicts source counts 
of ~1.7xl0 sr and ^1.8xl0 5 sr _1 for sources brighter than lOOmJy at 
150 urn and 180 urn, respectively. The number of sources found in this study 
clearly exceeds these predictions. 

Franceschini et al. || have presented similar models which include con- 
tributions from two galaxy populations: dust-enshrouded formation of early- 
type galaxies and late-type galaxies with enhanced star-formation at lower 
redshifts. The predicted source counts at 170pm are higher than in the model 
of Guiderdoni et al. j^] and thus in better agreement with our results. 

In Fig. 3 we show the predictions of these two models together with counts 
from other references. For this plot we have selected from Fig. || the values 
of the dotted histogram. 
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Fig. 2. Cumulative source counts at 90(i.m, 150 [im and 180(i.m. The dotted line 
gives all sources detected and the solid line sources detected at more than one wave- 
length. Areas were estimated according to the faintest source detected in a map. 
The dashed curves represent cases in which both the the source selection and the 
area determination were based on the local background fluctuations (see text). 
Fig. 3. Comparison with other ISOPHOT counts and models of galaxy evolution. 
Our source counts at 100 mjy and 200 mjy are shown together with the results of 
Puget et al. [^o| (circles), Kawara et al. Jl(]] (triangles), Oliver et al. [[l8| (diamonds) 
and Linden- V0rnle (stars). We present our results as horizontal lines that in- 
dicate the difference between the DIRBE calibration and the adopted ISOPHOT 
calibration (see text). The predictions of model E of Guiderdoni et al. fl] are shown 
with solid lines and the evolutionary model of Franceschini et al. Q with a dashed 
line. Dotted lines show predictions of no-evolution models (90u.m: Guiderdoni et al. 
Js]; 150u.m and 180u.m: Franceschini et al. |^|) 



